Field dependent quasiparticles in a strongly correlated local system II 
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We extend the renormalized quasiparticle description of the symmetric Anderson model in a mag- 
netic field H , developed in earlier work, to the non-symmetric model. The renormalized parameters 
are deduced from the low energy NRG fixed point for arbitrary field values. We find quasiparticle 
resonance widths, A a (H), which depend in general on the spin a as well as H. The low tempera- 
ture static properties can be expressed completely in terms of these parameters, which can also be 
used as inputs for a renormalized perturbation theory. We show that taking into account repeated 
quasiparticle scattering gives results for the longitudinal and transverse dynamic spin susceptibilities 
which are in very good agreement with those obtained from direct NRG calculations. 
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I, INTRODUCTION 

In an earlier paper— (hereafter referred to as I), we 
showed how the low energy behavior of the particle-hole 
symmetric Anderson impurity model^ can be described 
in terms of field-dependent quasiparticles. This model is 
characterized by the three independent parameters, Ed, 
the impurity level, A, the broadening of this level due 
to the hybridization with conduction electrons, and U, 
the interaction at the impurity site. In the absence of 
a magnetic field, it was shown earlier^^ that the low 
energy behavior can be described by an effective version 
of the same model with three corresponding renormal- 
ized parameters, A, and U. Subsequently, this ap- 
proach was extended to include a magnetic field H, and 
the parameters were then found to be field-dependent, 
£d,a(H), A a (H), and U(H). One way to calculate these 
parameters is from the low energy excitations of the nu- 
merical renormalization group (NRG) fixed point which 
was used in the earlier paper I, A strong magnetic field 
tends to freeze the spin fluctuations and leads to a de- 
renormalization of the quasiparticles. On increasing the 
field from zero in the strong coupling case the parameters 
for the quasiparticles slowly revert to their uncorrelated 
mean field values in the extreme high field limit. 

The renormalized parameters are not just a convenient 
way of describing the low energy behavior; they com- 
pletely specify the model. A renormalized perturbation 
theory (RPT) can be set up in which the free propagators 
correspond to fully dressed quasiparticles^. This for- 
malism is particularly effective for describing the Fermi 
liquid regime, as only diagrams up to second order have 
to be taken into account to obtain asymptotically exact 
results for the T = susceptibilities, and the leading T 2 
term in the conductivity (I). This perturbation expansion 
is not restricted to the low energy and low temperature 
regime, and can be used for calculations on all energy 
scales. We have shown that a very good description of 
the T = spin and charge dynamics for the Anderson 
model in the Kondo regime can be obtained by summing 
the RPT diagrams for repeated quasiparticle scattering-. 
The results give an accurate description of the spin and 



charge susceptibilities for arbitrary magnetic field val- 
ues H, and for frequencies w extending over a range sig- 
nificantly larger than the Kondo temperature Xk. The 
Korringa-Shiba relation^ and the sum rules for the spec- 
tral density are satisfied. 

In this paper we show that this approach can be ex- 
tended to the non-symmetric Anderson model. There are 
significant differences in this the parameters ac- 

quire a spin dependence, and formulae given earlier have 
to be generalized. 



II. THE NON-SYMMETRIC ANDERSON 
MODEL IN A MAGNETIC FIELD 

The Hamiltonian for the Anderson model^ is 



-Ham = e^adld^ + Un 



d,T"<U 



(1) 



k.fT k,a 

where Ed,a — £d — ag^Hjl is the energy of the local- 
ized level at an impurity site in a magnetic field H, U 
the interaction at this local site, and Vk,a the hybridiza- 
tion matrix element to a band of conduction electrons 
of spin a with energy Ek, a — c<? c /ZB-ff/2, where g c is the 
g-factor for the conduction electrons. When [7 = the 
local level broadens into a resonance, corresponding to 
a localized quasi-bound state, whose width depends on 
the quantity A CT (u) = nj^k \ v k,a\ 2 S(Lu - £fc, CT ). For the 
impurity model, where we are interested in universal fea- 
tures, it is usual to take a wide conduction band with a 
flat density of states so that A CT (ix>) becomes independent 
of u>, and can be taken as a constant A CT . In this wide 
band limit A CT (w) will be independent of the magnetic 
field on the conduction electrons, so we can effectively 
put g c = 0. When this is the case A CT is usually taken to 
be a constant A independent of a. 

In the renormalized perturbation theory approach&i 
we cast the corresponding Lagrangian for this model 
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FIG. 1: (Color online) Plots of the renormalized parame- 
ters, A T (»/A, A|(/i)/A, £ d (ft)/7rA, f/(/i)/f7, ^(ft), for the 
asymmetric Anderson model, with 7rA = 0.1, U/ttA = 2 and 
Ed/ it A = —0.3, as a function of the logarithm of the magnetic 
field h/nA. The ratio A|(/i)/A T (/i) is also shown. 



jCam(£(1,(t, A, J7) into the form, 

^AM^d.er, A, {/) = £AM(£d,cr, A CT , [/) + £ c t(Al, A 2 , A3), 

(2) 

where the renormalized parameters, £d, CT and A CT , are de- 
fined in terms of the self-energy S CT (o;) of the one-electron 
Green function for the impurity state, 



(3) 



W - Eda + «A - Scr(w) 

and are given by 

Ed,a =Z CT (£d,<r + S CT (Q)), A CT = Z CT A, (4) 

where z CT is given by z a = 1/(1 — S^(0)). The renormal- 
ized or quasiparticle interaction U, is defined in terms 
of the local total 4-vertex r-|-|(a;i, u>2, 0J3, 0J4) at zero fre- 
quency, 



[7 = ^1^(0,0,0,0). 



(5) 



It will be convenient to rewrite the spin dependent quasi- 
particle energies in the form, e^a = £d{h) — ahrj(h), 
where 



Sd 



(6) 



where id(h) and fj(h) are both even functions of the mag- 
netic field h = gp^H/2. 

The renormalized perturbation expansion is in pow- 
ers of the renormalized interaction U for the complete 
Lagrangian defined in equation ([2]). The counter term 
part of the Lagrangian £ c t(Ai, A2, A3) essentially takes 
care of any overcounting. The parameters, Sd.a, A CT and 
U, have been taken to be the fully renormalized ones. 



The most effective way of estimating them in this case 
is from the energy levels of a numerical renormalization 
group (NRG) calculation^. However, the approach it- 
self is independent of the NRG, and the parameters can 
be deduced in other ways, for instance from experiment. 
The counter term parameters, Ai, A2 and A3, are required 
to cancel any further renormalization, and are completely 
determined by this condition^. 

In figure Q] we display some typical results for renor- 
malized parameters as a function of the magnetic field 
on a log scale. For this plot we have taken bare pa- 
rameters Ed/VA = —0.3 and U/nA = 2, correspond- 
ing to a impurity occupation in the absence of a field, 
( n d,a) = n(0)/2 ~ 0.35. With these values the param- 
eters are strongly renormalized for h = 0. The over- 
all trend as a function of h is very similar to that for 
the particle-hole symmetric case (I) in the strong cou- 
pling regime, with the parameters converging to the mean 
field values in the limit h — » 00. We do see, how- 
ever, that the resonance widths become spin dependent, 
A|(/i) 7^ A;(/i), except asymptotically as h — ► and 
h ► 00. We note that, though Af(/i) increases mono- 
tonically with increase of h, A±(h) initially decreases. In 
this case, where the impurity level is less than half-filled, 
the ratio Af (/i)/Aj p (/i) > 1. This ratio is reversed, so 
Af (h)/ A.i (h) < 1, when the impurity level is more than 
half-filled. 



III. STATIC RESPONSE FUNCTIONS 

Once the renormalized parameters have been deter- 
mined, the static impurity response functions can be 
evaluated directly by substituting into the relevant ex- 
act formula. This is a more direct way of calculating 
these response functions than the usual NRG way, which 
involves a subtraction procedure to estimate the impu- 
rity componen t 10 ' 11 ' 12 . The formulae in the absence of 
a magnetic field^, and for the symmetric model in the 
presence of a field, were given earlier (I). Here, we give 
the generalizations for the non-symmetric model. The in- 
duced magnetization M{h) is given by M(h) — gp,Bm(h), 
where 



m(h) = i(n dT - n dl ) = — ^ ^ a 



tan 



£dcr(fr) 



, (7) 



which can be derived from the Friedel sum rule. The 
longitudinal susceptibility Xl{h) (in units of (s^b) 2 ) is 
given by 

Xi(h) = 0.25(p T (0, h) + p x (0, h) + U(h)fr(0, typ^O, h)), 

(8) 

where p a {u>,h) is the free quasiparticle density of states 
given by 



1 



A a (h) 



tt ( W - idAh)) 2 + A2(ft) 



(9) 
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FIG. 2: (Color online) The induced magnetization m(h) as 
a function of the logarithm of the magnetic field h for the 
asymmetric Anderson model with the same set of parameters 
as given in figure [T] The dashed curve is that calculated 
from the direct evaluation of the occupation values from the 
NRG ground state, and the full curve is that deduced from 
the renormalized parameters in equation (|7|). Also shown is 
the average occupation n(h)/2 as calculated from the NRG 
ground state (dashed curve) and the quasiparticle occupation 
values (full curve) as given in equation llllll . 



The corresponding transverse susceptibility Xt{h) (zero 
applied field limit in the transverse direction) is given by 



Xt(h) 



m{h) 
2h 



(10) 



The total occupation of the impurity site n(h) = (n^ + 
ndx) can be derived similarly, and is given by 



^) = l-5=E 



. tan" 1 (2^ 



(11) 



and the local charge susceptibility Xc(h) is given by 

Xc(h) = 0.25[p T (0, h) + ^(0, h) - U(h)fr(0, hjp^Q, h)}. 

(12) 

To illustrate the magnetic field dependence of these 
static response functions, we evaluate these formulae in 
a particular case using the field dependent renormalized 
parameters given in figure [TJ The results are compared 
with those obtained by a evaluation of static expectation 
values using the NRG. In all subsequent calculations in 
this paper, we concentrate on this same set of bare pa- 
rameters, Ed/fA = —0.3 and U/nA = 2. 

Shown in figure [2] is the result for m(h) deduced from 
equation ([7]), compared with results obtained by the di- 
rect evaluation of the d-site occupation values in the 
ground state as determined from the NRG. There is a 
small but systematic difference, of the order of 2%, be- 
tween the two sets of results. This difference could be 
due to the fact that we assume an infinite bandwidth 
in the derivation of these formulae, whereas we take a 
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FIG. 3: (Color online) The longitudinal and transverse impu- 
rity site susceptibilities, XiQ 1 ) and Xt(h), as a function of the 
logarithm of the magnetic field h. XiQ 1 ) ' s calculated from 
equation ((8]) and Xt(h) = m(h)/2h. 



finite band width 2D, with D = 1 in the NRG calcu- 
lations. The corresponding estimates of the average oc- 
cupation number n(h)/2 as a function of magnetic field 
h are shown in the same figure. For this quantity the 
two sets of results are almost indistinguishable. In the 
extreme large field limit the average occupation of the 
impurity level tends to unity, as the majority spin level 
gets pulled further and further below the Fermi level, 
and the impurity becomes completely polarized. In this 
regime the average renormalized level, id, approaches the 
mean field value, id = £d + 0.5Un(h), rather than the 
bare value £d, while the other renormalized quantities 
approach their bare values. 

The longitudinal and transverse spin susceptibilities, 
Xl{h) and Xt{h), are plotted in figure [3] as a function of 
the logarithm of the magnetic field h. Xl(h) ls calculated 
from equation J8]) and Xt{h) = m{h)/2h as calculated 
from equation |(7j). They should asymptotically converge 
to the same result in the limit h — > 0. There seems to be 
a very small discrepancy, of about 1%, between the two 
estimates in this limit. This could be due to evaluating 
the ratio in equation (flQ| for very small values of the 
field. 

We now look at the behavior of the corresponding dy- 
namic response functions. 



IV. DYNAMIC RESPONSE FUNCTIONS 

A. Single particle spectra 

It is of interest to see how well the free quasipar- 
ticle density of states given in equation ([9]) compares 
with the spectral density p a (co) calculated from the NRG 
for each spin type in the presence of a field. The 
NRG spectra are calculated from the self-energy follow- 
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FIG. 4: (Color online) The spectral density for the spin up 
(majority) electrons pf(cj,h) (full curve) and spin down (mi- 
nority) electrons pi(io,h) (dashed curve) at T = as a func- 
tion of w for h/ffA = 0.001 (upper panel) and H/tvA — 0.04 
(lower panel). The two curves on the upper panel are almost 
coincident on the scale shown. Also shown are the correspond- 
ing quantities derived from the free quasiparticle densities of 
states, Z]Pf(io) (dot-dashed) and zipi(u) (dotted). 



ing the prescription given ir>i£, and we have also used 
the improved metho d 14 ' 15 based on the complete Anders- 
Schiller basis^. 

In figured] (upper panel) we show p-f(w) and pj.(w), as 
calculated from the NRG, for the field h/irA = 0.001. 
We see that in this relatively weak field that there are 
only a small differences between the two spectral densi- 
ties; they are almost coincident on the scale shown. The 
corresponding quantities derived from the free quasipar- 
ticle densities of states, z^pt(u>) and zipi(u>), are shown 
in the same figure. 

In figured] (lower panel) we make a similar comparison 
for a large field h/irA = 0.04, and the same set of bare 
parameters. The polarization is much stronger for this 
higher magnetic field value and there is now a marked 
difference between the two spectral densities for the two 
spin types. The spectra derived from the free quasiparti- 
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FIG. 5: Repeated scattering of a quasiparticle with spin f 
and a quasihole with spin J. via the effective interaction Ut. 



cles, can be seen to describe the spectrum in the immedi- 
ate vicinity of the Fermi level for both spin types. They 
cannot describe the difference in the heights of the peaks 
of the two spin types, because in these simple Lorentzian 
formula the height of peaks in the quasiparticle weighted 
spectra are independent of H and the spin and given by 
1/nA. Note that the description of the low frequency 
range can be extended within the renormalized pertur- 
bation theory frame work by including a self-energy cor- 
rection in the quasiparticle density of states as shown in 
the particle hole symmetric case^. 

In the next section we look beyond the simple free 
quasiparticle picture and take account of the leading 
quasiparticle scattering terms. 



B. Dynamic spin susceptibilities 

The leading corrections to the two particle dynamic 
response functions to the free quasiparticle picture in 
the renormalized perturbation theory arise from the dia- 
grams with repeated quasiparticle scattering. In figured] 
we illustrate this type of diagram for the scattering of a 
free up spin quasiparticle with a free down spin quasihole. 

Such a diagram contributes to the transverse dynamic 
susceptibility Xt(v,h) in the presence of the magnetic 
field. We derived a quasiparticle interaction term U(h) 
earlier, but this is not the interaction term we need in 
considering this diagram, as this corresponds to the to- 
tal 4-vertex, (0, 0, 0, 0), which has implicitly included 
this scattering for u) — 0. What we use is the irreducible 
particle-hole interaction at zero frequency in this chan- 
nel, which we denote by Ut{h). The dynamic transverse 
susceptibility x*( w j h), is then given by 



l-U t {h)x n {uj,h) 



(13) 



where xu( w 7 h) is the transverse spin susceptibility of the 
free quasiparticle. The analytic expression for Xfj.(w,/i) 
is given in the appendix. 

To find the unknown quantity Ut(h) we use the fact 
that this series for ui = gives the static transverse sus- 
ceptibility which we have calculated already, and is given 
in equation (jTDJ) . Hence, comparing (fT3|) for u> = with 
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equation lfT0|) yields 



Ut(h) 



Ah 2 



(At-a x ) 2 



m{h) 



Ah 



(A T -A;) 
27rm(/i) 



In 



(ft-e d ) 2 +A^ 



(14) 

where h — hfj(h), and m(h) is given in equation 0- 
This expression simplifies in the case, Aj = A^, to give 
Ut{h) = (h — h)/m(h), which is the same as that used 
earlier—. 

Now we test how well these formulae describe the spin 
dynamic response over the relevant frequency range. To 
do this we compare the predictions based on the RPT 
formula lfT3|) with those obtained from a direct NRG eval- 
uation. 

The imaginary part of the transverse dynamic suscep- 
tibility is shown in figure [6] (upper panel) for a magnetic 
field value h/irA = 0.001. The dashed curve is that from 
a direct NRG evaluatio n 18 ' 19 and the full curve is that 
calculated using equation (fl3|l . with the corresponding 
renormalized parameters. There is very good agreement 
between the two sets of results. 

Results for the imaginary part of Xt(w, h) for the field 
value h/nA = 0.04 are shown in figure [6] (lower panel). 
In this stronger field one of the peaks is suppressed while 
the other peak is enhanced. The peak positions are in 
good agreement, and the slightly broader peak from the 
NRG data can be attributed to the logarithmic broaden- 
ing used in the direct NRG evaluation. There is a sum 
rule, that the total spectral weight is equal to —2m(h), 
which is satisfied precisely both in the RPT result and 
also in the NRG calculation, as we have used the im- 
proved prescription for the response functions based on 
the complete Anders-Schiller basis. 

In a similar way, we can calculate the longitudinal 
dynamic susceptibility Xl(u,h) from repeated scattering 
of the spin up and spin down quasiparticles. The ir- 
reducible zero frequency interaction in this channel we 
denote by Ui(h), and determine it from the condition, 
Xi(0, h) — xi(h), using equation J8]). The results for 
Xi(u,h) is 

Xi{u,h)= (15) 

XTt(^> h ) + Xu(^> h ) + 4t/i( fe )XTT(^> h )xa(^, h ) 
2(l~4U?(h)x n (u;,h)x ll (u;,h)) 

where the analytic expression for x<t<t(w, h) is given in the 
Appendix. Equation (fT5"l) is a generalization to our earlier 
result for the particle-hole symmetric model. The zero 
frequency irreducible particle-hole vertex Ui(h) in this 
scattering channel is given by 



Ui(h) = 



-l + j [i + u 2 (h)(pi ip(h)(pi -~ Pl ) 2 



(16) 

where we have simplified the notation, p a (0, h) = p a . In 
the absence of a magnetic field, or with a magnetic field 
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FIG. 6: (Color online) The imaginary part of the transverse 
dynamic susceptibility, xt(w, /i) at T = as a function of tu, 
for h/irA — 0.001 (upper panel) and h/nA — 0.04 (lower 
panel). The dashed curve is calculated from a direct NRG 
calculation and the full curve from equation Ill3ll with the 
renormalized parameters. 



for the particle-hole symmetric model, p| = pj,, the result 
simplifies to U t (h) = U(h)/[l + U(h)p(Q,h)], which is the 
value used in the earlier work^. 

The imaginary part of the longitudinal dynamic spin 
susceptibility in the weak field case, h/irA = 0.001, is 
similar to that for the transverse case shown in figure [6] 
(upper panel), apart from an overall factor of 2, so we 
do not show he results here. The stronger field case, 
h/nA = 0.04, is shown in figure 

Due to the stronger magnetic field, the heights of the 
peaks are slightly reduced as compared with the results 
for h/nA — 0.001, but the overall features are very sim- 
ilar. The Korringa-Shiba relation does not hold for the 
model without particle-hole symmetry in the presence of 
a magnetic field. Where it does hold, in the absence of 
a field, or with particle-hole symmetry, the renormalized 
perturbation expression satisfies it exactly 8 -. 
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FIG. 7: (Color online) The imaginary part of the longitudinal 
dynamic susceptibility, xifajh) at T = as a function of u 
for (h/nA = 0.04). 



V. SUMMARY AND DISCUSSION 



tral density p CT (w), except near the Fermi level. Hence 
in the strong correlation regime, when z a is very small, 
it would appear to be more appropriate to interpret the 
quasiparticle as describing a spinon excitation. Because 
the dominant low energy excitations are spinons, as is 
known from the Bethe ansatz solutions for the Kondo 
model 2l i»2i, then it is not so surprising that they provide 
a good description of the spin dynamics. 



Acknowledgement 

We wish to thank N. Dupuis, D.M. Edwards, W. Roller, 
D. Meyer and A. Oguri for helpful discussions and W. 
Roller and D. Meyer for their contributions to the devel- 
opment of the NRG programs. One of us (J.B.) thanks 
the Gottlieb Daimler and Rarl Benz Foundation, the Ger- 
man Academic exchange service (DAAD) and the EP- 
SRC for financial support. 



In this paper we have extended our earlier work I, 
where we described the low energy behavior of the sym- 
metric Anderson model in a magnetic field h in terms of 
field-dependent renormalized quasiparticles, to the non- 
symmetric Anderson model. The main new feature that 
emerges is the dependence of the quasiparticle peak res- 
onance width A a (h) on the spin type a as well as on 
the value of the magnetic field. The T — spin and 
charge susceptibilities can be expressed as exact formu- 
lae in terms of the parameters A a (h), id,a(h) and a lo- 
cal field dependent interaction between the quasiparticles 
U(h). Based on these parameters the low temperature 
behavior of the model, such as the susceptibilities and 
conductance, in the presence of the field could be calcu- 
lated as described for the symmetric model in I. It was 
also demonstrated earlier that an excellent description of 
the low energy spin dynamics can be obtained for the 
symmetric model- based on approximate formulae which 
take into account repeated quasiparticle scattering in the 
RPT. Here these results have been generalized to the non- 
symmetric model for the transverse and longitudinal spin 
susceptibilities, which again agree remarkably well with 
those obtained from a direct NRG calculation. 

What may seem surprising at first is that the results 
for the dynamic susceptibilities, based on the repeated 
quasiparticle scattering, agree so well over the full range 
of u, whereas the free quasiparticle density of states only 
describes the single electron spectral densities well only in 
the immediate region of the Fermi level. However, as dis- 
cussed in I, the quasiparticle density of states p a (ui), sur- 
vives in the limit z a — > 0, when it becomes a delta func- 
tion describing a free spin. For very small z a the prod- 
uct Zapa^ui), makes very little contribution to the spec- 



VI. APPENDIX 

The free quasiparticle dynamic susceptibility X<*,<j' ( w ) 
for the impurity model in the wide band limit, Aj = Aj,, 
were given earlier 8 . Here we give the more general results 
for A T ^ A ; , 



-1 A CT 



E Mi 



a=— 1,1 



aid.a- + iA a 



for li > 0, and for u> = 0, 

W(0) = p tr (0). 



(17) 



(18) 



The values for to < follow from the fact 
that Rex CTj(T (w) = Rex CTjer (— u) and lmx Uj(T {u) = 
-Imx CT , CT (-w). For a' a, 



i/2n 



U - id,] - i&] 
(w + e d ,i -£d,t +* A T -* A l) V -i^l-Cdd 



-In 



i/2tt 



-In 



g + £d,i - »A X 
(uj + id, i - e d ,] - iA T + iAi) y -iA]+id,-\ 
-i/2-K 

■ = = — x 

(w + erfj - id,] +«A T +iA;) 



In 1 



zA T - id, 



In 1 



* A 1 +?d,l 
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